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ABSTRACT 

The long-duration GRB050713a was observed by the MAGIC Telescope, 40 
seconds after the burst onset, and followed up for 37 minutes, until twilight. The 
observation, triggered by a SWIFT alert, covered energies above ~ 175 GeV. 
Using standard MAGIC analysis, no evidence for a gamma signal was found. As 
the redshift of the GRB was not measured directly, the flux upper limit, estimated 
by MAGIC, is still compatible with the assumption of an unbroken power-law 
spectrum extending from a few hundred keV to our energy range. 

Subject headings: gamma rays: bursts — gamma rays: observations 
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1. Introduction 

Observations of high-energy photons from gamma-ray bursts (GRBs) have much con- 
tributed to a deeper understanding of their nature. The 7-ray emission observed by the 
EGRET detector (Hurley et al. 1994) suggests a power-law spectrum extending up to GeV 
energies. This favours an optically thin emission region and a non-thermal origin of the 
bursts. As the excessive pair production could be suppressed in the presence of relativistic 
jets (Goodman 1986; Paczynski 1986), it was concluded that relativistic beaming could play 
an important role for GRBs (Meszaros & Rees 1993). However, other models also point 
towards the presence of a strong thermal component in the GRB spectra (Ryde 2004). 

The observation of 7-rays at highest energies is expected to have an important impact on 
the modelling of the emission processes, in particular of the early and late afterglow phases 
of GRBs. EGRET measurements generally showed the presence of a hard, long-duration 
component (Dingus 1995), consistent with a simple extrapolation of the MeV spectrum into 
the high-energy 7-ray regime. Recently, an additional, delayed high-energy component of 
GRB970417 was found with the TASC detector of EGRET (Gonzalez et al. 2003). Several 
models predict GeV-TeV emission lasting up to the early afterglow (Pe'er & Waxman 2004; 
Dermer & Atoyan 2004). Due to the extremely high energies attainable inside relativistic 
jets, GRBs are potential sources of Very High Energy (VHE) cosmic rays (Waxman 1995; 
Vietri 1995), that can produce in their turn hadronic showers containing VHE 7-rays. Other 
theoretical models predict no emission above few MeV (Lazzati et al. 2004) or predict strong 
emission up to GeV, but no emission above 10 GeV (Stern & Poutanen 2004). Therefore, 
measurements at this energy range can be used to test all these competing models. However, 
as most of the observed GRBs occur at large redshift, strong attenuation of the VHE 7 flux is 
expected, as a result of the interaction with low energy photons of the Metagalactic Radiation 
Field (MRF) (Nikishov 1961; de Jager & Stecker 2002). The knowledge of the redshift is, 
therefore, important for a precise interpretation. On the other hand, a detection of VHE 
7-rays provides an indirect — and model dependent — upper limit of its redshift, if some 
knowledge of the MRF is assumed. 

Several observations of GRBs at energies above hundred GeV have been attempted 
(Goetting et al. 2003: GCN #1007 ; Zhou et al. 2003), without showing any indication of 
a signal. This is due to a relatively low sensitivity, as in satellite-borne detectors, and/or 
a high energy threshold, as in previous generation of Cherenkov telescopes or in particle 
detector arrays. Up to now, only upper limits on the prompt or delayed emission of GRBs 
were set by Whipple (Connaughton et al. 1997), MILAGRO (see Atkins et al. 2005, and 
refs. therein), and STACEE (Jarvis et al. 2005). STACEE, in the same energy region as 
attainable by MAGIC, was able to follow GRB050607 after 3'11" for 1150 s and set an 



-4- 



upper limit of its flux as $(> fOO GeV) < 4.1 x 10~ 9 cm^s" 1 « 6 C.U. (Cra& Ms: 
1.5 x 10~ 6 x £(GeV)~ 2 - 58 ph • cm' 2 s^CcV" 1 ). 

The situation may change with the new generation of Cherenkov telescopes, which 
achieve a better flux sensitivity and a lower energy threshold. Nevertheless, as their small 
field of view does allow prompt observations only by serendipitous detection, they have 
to rely on an external triggering, as the one provided by the automated satellite link to 
the Gamma-Ray Burst Coordinates Network (GCN) 1 , which broadcasts the coordinates of 
events triggered and selected by dedicated satellite detectors. 

Among the new Cherenkov telescopes, MAGIC (Mirzoyan et al. 2005) is best suited for 
the detection of the prompt emission of GRBs, due to its low energy threshold, large effective 
area, and, in particular, its capability for fast slewing (Bretz et al. 2003). The low trigger 
threshold, currently 50 GeV at zenith, should allow the observation of GRBs even at large 
redshift, as lower energy radiation can effectively reach the Earth without interacting much 
with the MRF. Moreover, in its fast slewing mode, MAGIC can be repositioned in < 30 s 
to any position on the sky: in case of a Target of Opportunity alert by GCN, an automated 
procedure takes only few seconds to terminate any pending observation, validate the incom- 
ing signal and start slewing toward the GRB position. Extrapolating BATSE observed GRB 
spectra to VHE with an unbroken power-law of power index from the BATSE catalogue, 
MAGIC is predicted to detect about one GRB per year at a 5<r level (Galante et al. 2003). 

In this letter, we report on the analysis of data collected on GRB050713a during its 
prompt emission phase and for the following 37 minutes. 

2. MAGIC observation 

On 2005 July 13 at 4:29:02 UT the BAT instrument on board SWIFT detected a burst 
located at RA 21 h 22 m 09*53 DEC +77°04'29 , .'50±3' (Falcone et al. 2005: GCN #3581). The 
MAGIC alert system received and validated the alert 12 s after the burst, and data taking 
started 40 s after the burst original time (T ) (Galante et al. 2005: GCN #3747). 

The burst was classified by SWIFT as a bright burst with a duration of T 90 = 70 ± 10 s. 
The brightest part of the keV emission occured within T + 20 s, three smaller peaks followed 
at T + 50 s, T + 65 s and T + 105 s, while a pre-burst peak took place at T — 60 s. (see 
figure 1). The spectrum, over the interval from T — 70 s to T + 121 s, can be fitted with a 
power-law with photon index —1.58 ±0.07 and yields a fluence of 9.1 x 10 _6 erg • cm -2 in the 
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15^-350 keV range (Palmer et al. 2005: GCN #3597). The burst triggered also Konus-Wind 
(Golenetskii et al. 2005: GCN #3619), which measured the spectrum of the burst during the 
first 16 s, that is the duration of the first big peak as reported by SWIFT. 

2.1. Data set and analysis 

In the local coordinate system of MAGIC, GRB050713a was located at an azimuth angle 
of —6° (near North) and a zenith angle of 50°. The sky region of the burst was observed 
for 37 min, until twilight (ON data). Between T + 665 s and To + 686 s, data taking was 
interrupted for technical reasons. A total amount of 258250 atmospheric showers, mainly 
background, were recorded. In order to evaluate the contamination of background in the 
data, the GRB position was observed again 48 hours later (the so-called OFF data). 

Data were analyzed using the MAGIC standard software (Bretz et al. 2005; Gaug et al. 2005). 
For optimizing 7/hadron separation, we simulated 10 5 7 events with zenith angle ranging 
between 47° and 52°, energy greater than 10 GeV and an energy distribution following a 
power-law spectrum of index /3 7 = —2.6. This sample was analyzed in the same way as the 
data, and was used for the calculation of the collection area, the sensitivity, and for the en- 
ergy estimation. After applying all selection criteria, the sample peaked at around 250 GeV, 
which we define as our telescope threshold at this zenith angle. 

Data are processed using the standard Hillas analysis (Hillas 1985; Fegan et al. 1997). 
Gamma/hadron separation is performed by means of Random Forest (RF)(Breiman 2001), 
a classification method that combines several parameters describing the shape of the image 
into a new parameter called hadronness, the final 7/hadron discriminator in our analysis. 
The simulated sample was used to optimise, as a function of energy, the cuts in hadronness. 
Also the energy of the 7 was estimated using a RF approach, yielding a resolution of « 30% 
at 200 GeV. 

The parameter alpha of the Hillas analysis, related with the direction of the incoming 
shower, is used to evaluate the significance of a signal. If the telescope is directed at a point- 
like 7 source, as a GRB is expected to be, the alpha distribution of collected photons should 
peak at 0°, while it is uniform for isotropic background showers. According to simulations 
the 7-signal at low energies may spread in a region defined conservatively by alpha < 30°. 
Figure 2 shows the alpha distributions for the GRB050713a and for OFF data, divided into 
three subsets of time covering 90 s, 5 min and 30 min, respectively. No evidence of excess in 
the signal region is seen. 
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2.2. Time analysis 

A second analysis searching for short time variable 7-ray signals from GRB050713a has 
been performed in the range 175 GeV < E < 225 GeV. Figure 3 shows the number of 
excess events during the first 37 minutes after the burst, in intervals of 20 s. The number of 
expected background events in the signal region (open circles), estimated from the number of 
events in the region with alpha > 30°, is constant indicating stable experimental condition. 
The number of excess events is stable and compatible with statistical fluctuations of the 
background. The same analysis was applied to the OFF data, with similar results. 

2.3. Flux Upper Limits 

Analysing the data collected during the prompt emission of GRB050713a between T + 
40 s and T + 130 s, we can set upper limits on its flux at 95% confidence level (see details 
in: Rolke, Lopez, & Conrad 2005). 

The upper limit can be used to constrain the prompt emission of the GRB in the VHE 
range. Since the observed spectrum is the convolution of the intrinsic spectrum and the 
MRF absorption, the limits on the former are thus necessarily model dependent. 

First of all, we assumed the GRB spectrum extends to GeV energies following the Band 
function (Band et al. 1993): after the energy break, estimated by Konus-Wind to be at 
~ 355 keV, the flux follows a power-law of spectral index (3 = —2.5, the mean value of the 
BATSE distribution (see Preece et al. 2000). In this hypothesis, we calculated the upper 
limit on the average flux in our energy range during the entire 90 s interval. These values 
are summarised in table 1, and the lowest two energy bins are shown in figure 4, together 
with the spectrum measured at lower energies by SWIFT and Konus-Wind. 

It has to be noted, however, that according to BAT data, only 10% 15% of the 
total burst fluence in the 100 keV region was released in the time window of the MAGIC 
observations. This fraction of the flux is plotted in figure 4 using a dashed line. Adopting a 
semi-empirical model for the cosmologically evolving MRF (Kneiske et al. 2004), we derived 
unfolded flux upper limits for various redshift values as shown in table 1. 

3. Conclusions 

MAGIC was able to observe part of the prompt emission phase of a GRB as a response 
to the alert system provided by the SWIFT satellite. No excess above 175 GeV was detected 
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neither during the prompt emission phase nor during the following 37 minutes. We derived an 
upper limit to the 7-ray flux between 175 and 1000 GeV. The observation window covered 
by MAGIC did not contain the first prominent peak detected at keV energies where the 
SWIFT and Konus-Wind spectra were taken. Upper limits are compatible with any naive 
extensions of the power-law spectrum up to hundreds of GeV. 

For the first time a Cherenkov telescope is now able to perform direct observations of 
the prompt emission phase of GRBs. Although strong absorption of the high-energy 7-ray 
flux by the MRF is expected at high redshifts, given its sensitivity to low fluxes and its fast 
slewing capabilities, the MAGIC telescope is expected to detect about one GRB per year, if 
the GRB spectra extend to the hundreds of GeV energy domain. 
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Fig. 1. — MAGIC excess event rate compared with SWIFT-BAT observations. The vertical 
line indicates the start of observations with the MAGIC telescope: the prominent peak seen 
by SWIFT-BAT occurred before MAGIC observations started. 
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Fig. 2. — Alpha distributions of events with 175 GeV < E < 225 GeV for three different 
time intervals starting at T = To + 40 s: 30 min (top), 5 min (middle), and 90 s (bottom). 
Dots refer to ON data, the line to OFF data. The vertical line bounds the region where we 
expect the 7 signal. 
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Fig. 3. — Full circles: number of excess events for 20 s intervals, in the 37 min window after 
the burst. Open circles: number of background events in the signal region. 
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Fig. 4. — Upper limits set by MAGIC on GRB050713a with no redshift correction applied 
(see text). The solid line is the flux measured by SWIFT averaged over the burst T 90 , and the 
energy break is estimated using Konus- Wind data. The dashed line represents the fraction 
of the flux emitted between T + 40 s and T + 130 s. 
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Energy 


Excess evts. 


Eff. Area 


Flux upper limit x 10 9 [erg cm 


2 s- 1 } 


[GeV] 


upper limit 


xlO 8 [cm 2 ] 


z = 


z = 0.2 


z = 0.6 


z = l 


175 - 225 


8.5 


1.7 


0.83 (7.6 C.U.) 


1.16 


3.42 


10.49 


225 - 300 


10.4 


3.4 


0.45 (4.8 C.U.) 


1.07 


4.63 


19.32 


300 - 400 


6.0 


5.3 


0.37 (3.8 C.U.) 


1.35 


13.20 


95.45 


400 - 1000 


4.3 


6.5 


0.13 (3.3 C.U.) 


0.68 


25.11 


293.18 



Table 1: MAGIC upper limit (95% CL) on GRB050713a between T + 40 s and T + 130 s 
(see text). Limits include a systematic uncertainty of 30% and have been corrected for the 
photon absorption by the EBL for different redshift values z. For z = (no correction 
applied) the flux upper limits are expressed also in Crab Units (C.U.). 



